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Abstract: The cross section for central diffractive Higgs production is calculated, for the LHC range 
I/-) of energies. The graphs for the possible mechanisms for Higgs production, through pomeron fusion and 
photon fusions are calculated for all possibilities allowed by the standard model. The cross section for 
central diffractive Higgs production through pomeron fusion, must be multiplied by a factor for the survival 
probability, to isolate the Higgs signal and reduce the background. Due to the small value of the survival 
probability ( 4 x 10~ 3 ) , the cross sections for central diffractive Higgs production, in the two cases for 
q pomeron fusion and photon fusion, are competitive. 
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1. Introduction 

The most promising process for observation of the Higgs boson at the LHC is central diffractive production 
of the Higgs, with large rapidity gaps ( LRG ) between the Higgs and the two emerging protons, after 
scattering. Namely, 



p + p -> P + [LRG] + H + [LRG] +p (1.1) 

The large rapidity gaps either side of the Higgs reduces the background, so that the Higgs signal will be 
easier to isolate in central exclusive production. Hence, this process gives the best experimental signature 
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for detecting the Higgs at the LHC, and is very interesting for experiments in the search for the Higgs boson. 



In this paper, two mechanisms are compared for central exclusive Higgs production; (1) 77 fusion, 
namely pp — > 77 — > H and (2) pomeron exchange, namely pp — > IPIP — > H, where P denotes a 
pomeron. For 7 7 fusion, there is no hard re-scattering of the photons to fill up the rapidity gaps, so that 
large rapidity gaps are automatically present. The motivation for considering IP IP — > H, is the observation 
that for gluon gluon fusion, namely the process gg — ► H the colour flow induces many secondary parton 
showers which fill up the rapidity gaps. Instead the process IP IP — > H is a colour singlet exchange, where 
the colour flow is screened, and the large rapidity gaps are preserved. 



However in the case of IP IP fusion, there are hard re-scattering corrections, giving additional inelastic 
scattering, which will give emission filling up the rapidity gaps (see ref. [0). To guarantee the presence of 
large rapidity gaps after scattering in central exclusive production, in the case of IPIP fusion, one has to 
multiply by the survival probability. This is the probability that large rapidity gaps, between the Higgs 
boson and the emerging protons, will be present after scattering. 



The motive of this paper, was driven by the potentially small value for the survival probability, namely 
< I S 2 \ >= 4 x 10~ 3 , calculated in ref. M for central exclusive Higgs production via IP IP fusion. This 
is an order of magnitude less than previous estimates, namely < | S 2 \ > = 0.02 for LHC energies (see ref. 
[Q]). Hence, previous calculations in ref. [Q] of the cross section for central exclusive Higgs production via 
IPIP fusion, which included a factor for the survival probability of < | S 2 \ > = 0.02, gave for the exclusive 
cross section crf^ (pp — > p + H + P) = 3fb. Since the survival probability is predicted in ref. Q to be 
an order of magnitude less, it follows that the cross section crf^p (pp — ► p + H + P) for central exclusive 
Higgs production will also be an order of magnitude smaller. It also follows that this cross section will be 
competitive with the cross section for central exclusive Higgs production via 7 7 fusion, which is predicted 
in ref. f|] to be a^ c (p + p^ p + H + p) = O.lfb. 



To illustrate that for the process of Eq. ( |1.1| ), the cross sections will be competitive for 77 fusion 
and IP iPfusion, it is instructive to consider the diagrams for these processes. The notation used for the 
couplings in the standard model are the following. 
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coupling constant 


expression 


value ( GeV / c 2 ) 


C^em 


£_ 

An 


~ 137 


G F 




1.17 x 10~ 5 


m w 


\vg w 


80 


M H 


V2Xv 2 


120 


a s {Ml) 


9 2 s(M 2 H ) 
An 


0.12 



where the mass of the Higgs boson, is derived from the vacuum expectation value of the Higgs SU (2) 



weak isodoublet, which is 



vev 




, and and A are parametres in the 



Higgs potential, which is introduced into the standard model in the spontaneous symmetry breaking of 
SUl ( 2) x Uy ( 1) — ► Uem ( 1) ! which is responsible for giving rise to the W and Z boson masses. In the 
case for 77 fusion shown in Fig. ||, there are four vertices proportional to a em coupling the photons to the 
two protons and the photons either side of the subprocess for 77 — > H. A factor of g 2 ^ = 4\/2Gf vr? w 
should also be included, to account for the weak coupling of the Higgs to the sub-process, depicted by the 
shaded area in Fig. 0. So it is expected that the cross section will be proportional to 



<^ C (p + P-> p + H + P) oc 4V2G F m 2 w ai m = 0.6fb (1.2) 

where units are defined as 1 GeV~ 2 = 0.3893 mb . In central exclusive Higgs production for the case 
for 77 fusion, all the couplings of the photons shown in Fig. || are known constants, namely they are 
proportional to a em . Hence, the cross section for this diagram can be calculated exactly. On the other 
hand, when considering central exclusive Higgs production for the case for JPJP fusion shown in Fig. |l], 
the couplings of the gluons are not constants. In Fig. |l] there are four gluon couplings with the protons, 
giving a contribution to the cross section proportional to ( Q 2 ) , where Q 2 is the momentum transferred 
along the pomeron, and it is assumed that a s (Q 2 ) ~ 0.2. There are also two gluon couplings with the 
subprocess for JPJP — > H, giving a contribution to the cross section proportional to a 2 ( Mf^j . Taking 
the mass of the Higgs to be Mh ~ lOOGe^, then it is expected that a s (Af^-) ~ 0.12. Also in the case 
of JPJP fusion, a factor of g 2 ^ = 4^/2Gf should also be included, to account for the weak coupling 
of the Higgs to the sub-process shown in Fig. |l|. For JPJP fusion, this subprocess is the quark triangle 
subprocess shown in Fig. ^. Since the gluon itself couples weakly to the Higgs boson, only the contribution 
of the quark triangle subprocess is taken into account. The amplitude for the quark triangle subprocess, is 
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derived in section ^1 for the electromagnetic case, where a s replaces a em . After multiplying by a factor 
for the survival probability, which includes the survival probability < | 5 , ^ arc j| > = 4 x 1CP 3 which takes 
into account hard re-scattering of the pomeron, and < | <Sg fJ > = 5 x 10~ 2 which takes into account soft 
re-scattering of the pomeron, then one obtains for the exclusive cross section for central exclusive Higgs 
production, in the case of PP fusion the value 



a e ^(p + p^ p + H + p) oc 4V2G F m 2 w a 2 s (M 2 H ) a 4 s (Q 2 ) < |5j ard | ><|S s 2 oft | >= 0.9fb (1.3) 

Comparing the estimates of Eq. ( |1.2| ) and Eq. ( |1.3| ) , it is expected that (p + P — ► P + H + p) and 
aprp(p + p— * p + H + p) will be competitive. This is the motivation for this paper which re-examines 
(p + p — > p + H + p) for electromagnetic Higgs production. 

This paper is organised in the following way. In section [2L the details of the calculation of the cross 
section for central exclusive Higgs production in the case for PP fusion is explained in detail. The cross 
section which is obtained, is multiplied by the factor for the survival probability in ref. [jlj], to give the 
exclusive cross section. In section ||, the cross section for central exclusive Higgs production, in the case 
of 77 fusion is calculated. The mechanism 77 — ► H proceeds via the fermion triangle and Boson loop 
sub-processes illustrated in Fig. ||| and Fig. [|. The total cross section for central exclusive Higgs production, 
in the case of 77 fusion is calculated by taking the sum over all the contributions for these sub-processes 
for the 77 — > H mechanism. Finally, in the conclusion the results found in section [2] and section ||, are 
compared. 

All calculations in this paper, are based on the Feynman rules for the standard electro- weak theory, 
given in Fig. 0, which is to be found in the appendix. 

2. Double diffractive Higgs production at the LHC 

In this section, an explicit expression is derived for the Born amplitude for double diffractive Higgs pro- 
duction (see Fig. |l|). The Higgs couples weakly to the gluon, so the main contribution comes from the 
quark triangle subprocess (see Fig. ||), and all six flavours of quarks are taken into account. The second 
t channel gluon in Fig. [l] is included. This is because of the large rapidity (LRG) gap between the Higgs 
and the proton, which demands that in the t channel, one has colorless exchange. Indeed, if the LRG 
wasn't present between the protons, then the Higgs could simply be produced from gluon gluon fusion in 
a single channel. However, the colour flow induced by a single channel exchange process, could produce 
many secondary particles. These secondary particles could fill up the LRG. To screen the colour flow, it 
is necessary to exchange a second t channel gluon. At lowest order in a s , this gluon couples only to the 
incoming quark lines. The Born amplitude for double diffractive Higgs production by gluon exchange, is 
given by the expression || ||] 
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Figure 1: Double diffractive Higgs production in the Born approximation 



4 2s 



9M 2 



Aki-k 



d 2 Q fcir • k 2T 
k-^ 



8a 2 s (Q 2 



9M 2 



A- 



M 



H 



d 2 Q krr ; far 
(1^^ k-^ k>2 



8a 2 s (Q 2 ) 



(2.1) 



where k\- k 2 = has been used. In Eq. (|2.1| ), the Weizsacker - Williams approach, explained in ref. 
H has been used. The factor A k\ ■ k 2 is the amplitude for the quark triangle subprocess of Fig. |2|, where 
A takes the value, || 0, ||] 



.4 



7T 



(2.2) 



The Born amplitude shown in Fig. [I], is extended to proton scattering instead of just quark scattering 
here. The typical momentum transferred, t% = ( q\ — Q) 2 and t 2 = (q 2 — Q) 2 , are rather small and of the 
order of t, where b is the slope of the gluon - proton form factor, and can be estimated to be b = 5.5GeV 2 . 
Therefore, one takes into account the proton couplings to the gluon ladder (see Fig. |l|), by including the 
proton form factors in the gaussian form 



cxp 



1, 

— 6*i 
2 



-bt 2 
2 



where h = (qi — Q) and t 2 = ( q 2 - Q)' 



(2.3) 
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Higgs 



q denotes the sum over all flavour of quarks 



Figure 2: Quark triangle subprocess for Higgs production in gluon gluon fusion 

to describe the dependence on the transferred momentum. If the momentum transfer is small, it can 
be assumed that k\ ~ £?2 ~ Q- Therefore, with the definition of Eq. ( |2.3| ), t\ — ► and t2 — > 0. Hence, 
the proton form factors e~2 btl and e~z bt2 tend to 1, and can be taken outside the integral, and the Born 
amplitude behaves as 

2 -Ae-^ e -^ J ^?8a s (Q 2 ) (2.4) 

To consider the exclusive process only, with the condition of the LRG, bremsstrahlung gluons must be 
suppressed. The bremsstrahlung gluons are shown by the dashed lines in Fig. [l], which are suppressed by 
multiplying by the Sudakov form factor 

Fs = e ~S(kl,El) (25) 

F s is the probability not to emit bremsstrahlung gluons. S is the mean multiplicity of Bremsstrahlung 
gluons given as 

Secondly, evolution of BFKL ladder gluons between the two channels (see Fig. [I]), must be taken into 
account. For proton scattering instead of quark scattering, the naive coupling of the gluons to the external 
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quarks must be replaced by a coupling of the gluons to the external proton lines. To include both of these 
modifications, the naive gluon density for quarks is replaced by the density for protons by the following 
substitution 



^ -/(*.<« (2.7) 

where / (x, A; 2 ) is the un-integrated gluon density of the proton. After including the Sudakov form 
factor of Eq. (fO|), and the gluon density function of Eq. ( |2.7D , the amplitude in Eq. ( |2.1[) becomes 



Mjpjp (p+p - p+H+p) = Anh j ^L e - s ( k ^)f(x u Q 2 )f(x 2 ,Q 2 ) (2.8) 

f 2\ ^(71,2) lnl ) 

where for the gluon densities f { xxp) =2(<5) 7l ' 2 e \ X1 ' 2 J (2-9) 

where 71,2 are the anomalous dimensions and the numerical coefficient 2 in Eq. ( |2.9[ ) can be taken 



from MRST-2002-NLO parameterizations. Using Eq. (2.9), the integration over Q± and over 71 and 72 



can be evaluated. It turns out that the integrand of Eq. (|2.8|) has a saddle point given by In Q\ = 
In + ( 71 + 72 — 1) , and the essential values of 71 and 72 are close to |. Hence, the typical Q± is 
rather large, and depends on the mass of the Higgs. After integrating over 71 and 72 and q±, the final 
result for the amplitude Mjpjp ( p+p — > p+H+p) is derived in the appendix (see section A-l| ), and the final 



result is given in Eq. ( A-l-lq ) as 



in ^ 



VW ( fe+ln ^ +w » (i)l n 



SI 



Mjpjp (p+p -> p+H+p) =-2Att 4 s(-^) A- - 2 r — — (2.10) 



Now that the amplitude is known, ajpjp( p+p — > p+H+p) for central Higgs production can be cal- 
culated for the case of PIP fusion. To derive the cross section for exclusive central Higgs production, 
one has to multiply by a factor which takes into account the survival probability for large rapidity gaps 
< |>S'^ ar( j| > = 0.004, to suppress hard re-scattering. Therefore, the cross section <r|^p (p+p — > p+H+p) , 
for central exclusive Higgs production, without any further hard re-scattering for the case of JPJP fusion, 
takes the value 



< l^hardl > *w(p+P-> P+H+p) = ^ ( p+p ^ p+H+p) = 0.47fb (2.11) 
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3. Electromagnetic Higgs production at the LHC 



In the case of central exclusive Higgs production for the case of 77 fusion, shown in Fig. ^, there is no hard 
re-scattering to take into account, and all the couplings are known precisely. The shaded area in Fig. || 
depicts the subprocess for the mechanism 77 — ► H. The possible mechanisms are illustrated in Fig. || and 
Fig. H and their contributions to the amplitude are calculated in this section. Gauge invariance requires 
that the contribution of a subprocess for the mechanism 77 — * H, takes the form 



► photon 



Higgs 



q 2 'fv photon 



Figure 3: Diffractive Higgs production in single channel photon exchange. 



AT = A ( - ¥kgT ) (3.1) 



where A is a constant, depending on the particular subprocess. It turns out that for the case of when 
the subprocess is the fermion triangle shown in Fig. |], summed over all six flavours of quarks, and all three 
lepton flavours, then the expression for the amplitude takes the form of Eq. ( |3.1| ). However in the case of 
when the subprocess is one of the Boson loops shown in Fig. ^, the expression for the amplitude is not of 
the form of Eq. ( |3.ip . The correct statement is that the sum over all the amplitudes for the sub-processes 
shown in Fig. ||, gives a gauge invariant expression of the form of Eq. (|3.l|). 

The amplitude for the diagram of Fig. ^, where the process 77 — ► H proceeds via the sum over the 
fermion triangle shown in Fig. |], and Boson loops shown in Fig. |B|, is given by 

M 77 ( p+p - p+H+p) = - 4p£tf ( A{ v + A h \ (3.2) 

where A^ v = Af ^1^2 — ~^ L 9 flu ^j denotes the amplitude of the quark/anti-quark triangle subpro- 
cess, summed over all six quark flavours/anti-quark flavours (q,/q = u/u,d/d, s/s,c/c,t/t,b/b) and all 
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three lepton/anti-lepton flavours (L^ = e^ 1 , , ) , shown in Fig. |3|, and vl^ = ( q^q^ ■f-Q*"' ) 

denotes the amplitude of the sum over all the Boson loop sub-processes shown in Fig. |5[ At this point, the 
Weizsacker - Williams formula is used, as explained in refs. [5|, ||]. In this approach, the substitution 
PiP^A/tv = —jjTQi±q2±AiiiJ is used. In the notation used in this paper, q^ ± and q% ± denotes two dimen- 
sional vectors, in the plane transverse to the direction of the momenta of the two incoming protons p^ and 
P2- Hence, the amplitude of Eq. ( |3.2j ) can be written as 



M 77 (p+p- p+H+p) =-*^2L.4<£ L & L (A^ + A'J (3.3) 

<?1<?2 1V1 H 

In order to calculate the cross section, one has to integrate the squared amplitude over all the transverse 
momenta qi± and q-iL- For central exclusive Higgs production in the case of 77 fusion, it is required that 

the lower limits of integration are (</]j m ) 2 , (?2± n ) 2 = rn p\f^^'i where m p is the proton mass, which is 
assumed in this paper to be 1 GeV. The upper limits of integration, are taken from the electromagnetic 
form factors for the proton, namely G p {q 2 ) = —? — 1 2 \ ; from which the upper limits of the integration 

are derived to be {q^ x f = {qfl x f = 0.72. 



a- (p+p- p+ll + p) ^InJLiAf + Atftf I I CU) 



where y = In is the rapidity gap between the two incoming protons in Fig. ^, and m is the proton 
mass, assumed to be 1 GeV. This calculation is for central exclusive Higgs production at the LHC, where 
it is expected that = 14000 GeV, which gives for the value of the rapidity gap y = 19. 



3.1 The fermion triangle subprocess for Higgs production in 77 fusion 

Central exclusive Higgs production for the case of 77 fusion, can proceed through the subprocess 77 — ► 
fermion triangle — > H shown in Fig. |], where the fermions include the six flavours of quarks and anti 
quarks (u, d, s, c, t, b) and the three lepton and anti lepton flavours (e, fi, r). A derivation of the 
amplitude for the fermion triangle, can be found in refs. [^, ^, [l0| for the case where the mass of the 
fermion in the triangle is much larger than the Higgs mass. In this section, the amplitude of the fermion 
triangle is derived by taking the sum over all the fermions which could contribute, including the six quark 
flavours and the three lepton flavours. The H — > / / vertex coupling, is proportional to the mass mj of 
the fermion at the vertex, so that the sub-process amplitude of the fermion triangle will be proportional 
to the mass of the fermion in the triangle. The fermion masses are assumed to take the following values 
listed in the table below. 
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fermion 


mass ( GeV / c 2 ) 


quarks 




u 


d X 11) 


d 


D X 11) 


s 


1 Q 


c 


0.1 


t 


175 


b 


4.3 


fermions 




e 


5.11 x 1(T 4 


A» 


0.106 


r 


1.7771 



Hence, these values indicate that the most significant contribution to the amplitude will come from 
the top quark triangle. 




a; b) 



f~ = quarks q~ (u~ , d~ , s" , c , t~ , b ) and leptons L~ (e~fj~,T~) 
f + = anti -quarks q + (u + ,d + ,s + , c + , t + , b + ) and anti - leptons L + (e + ,/j + ,T + ) 



Figure 4: Fermion triangle subprocess for Higgs production through 77 fusion 
To derive explicitly the amplitude of the subprocess Fig. |], the labeling of momenta shown in the 
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diagram is used, and the following notation is introduced. 



D 1 =l 2 -m) D 2 = (I - q x f - m} D 3 = (I + q 2 ) 2 - m) (3.1.1) 

where rrif denotes the mass of the fermion which forms the triangle, which could be one of the quark 
flavours or one of the lepton flavours. Then the amplitude for the fermion triangle subprocess, summed 
over all possibilities of quark and anti-quark flavours, and lepton and anti- lepton flavours takes the form 



A 7 = £ 4,a em (V2G F ) 'mf —-A where E=E + E 



where Y^ q denotes the sum over all six quark flavours q = u , d , s , c , t , b and ^2 L = e „ T denotes 
the sum over all three lepton flavours L = e, /i, t. On the RHS of Eq. ( |3.1.2 ), d is the space-time 



dimension, and at the end f the calculation, cZ — ^ 4 is imposed. The reason for not specifying this in the 
beginning, is because it will be necessary to use dimensional regularisation to cancel divergences, which 
requires integration over d + e dimensions, in the limit that d — > 4 and e — > 0. Using the notation shown 
in the diagram of Fig. |] for the flow of momenta in the quark and anti - quark triangles, the trace term 
If U is given by 

If = Tr^ih + ql + mf)^ {l (T + m n v l T {lr-ql + rn f ))) 

+Tr (7* {-l x +ql + m f ) 7 V {-l a + mff-f {-l T - g? + m/))) 
= 8m f + 4ZT + 2 (l»t& - l L ' q%) - {q\ ■ q 2 + I 2 - m})) (3.1.3) 

where q\ and q 2 denotes four dimensional vectors. The first line on the RHS of Eq. ( |3.1.3| ) corresponds 
to the contribution given by the triangle formed by the fermions, ( i.e. quarks q = (u , d , s , c, t , b) 
and negatively charged leptons L~ = ( e~ , yT , t~) , and the second line corresponds to the contribution 
given by the triangle formed by the anti - fermions ( i.e. anti - quarks q = {u , d , s , c, t , b) and posi- 
tively charged leptons L + = ( e + , /i + , r + ) ). Introducing Feynman parametres to rewrite the quotient 
( D\ D 2 D3) 1 in a more convenient form, Eq. ( |3.1.2| ) simplifies to 

i l d d i r 1 , r 1 -' trnflf 



/ 

where A f = mj — x y 



(2tt) Jo Jo ' ([2_ A/ (x,y) 



and = -xq\ +yq 2 i (3.1.4) 

Note in the last step, it was assumed that q\ <C q\ Mfj^m 2 so that these terms can be ignored. 
From the kinematics shown in the diagram of Fig. [| and Fig. ||, it is clear that 2q\ ■ q 2 = Mg. The trace 
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term was given in Eq. ( |3.1.3 ) in terms of the unknown momentum in the fermion triangle in Fig. 
In terms of the new variable I 1 *, the trace term takes the form 



j-fJ,U 

.f 



m f 



Q1Q2 + ^ - y-sT[ qi ■ h (1 - 2x y) 



m 2 f + P 



(3.1.5) 



The details of the integration over the momentum I on the RHS of Eq. (3.1.4) are given in section A-3j 
of the appendix, (see Eq. ( [A-3-1 ) - Eq. ( A-3-4] )). Here, dimensional regularisation is used, a technique 



where one integrates over d + e dimensions, and afterwards d — > 4 and e — ► 
invariant terms, in the numerator of the integrand on the RHS of Eq. ( 3.1.4| ). 



following gauge invariant result for the RHS of Eq. ( 3.1.4 ). 



0. This removes non gauge 
In this way, one obtains the 



Af 



2a em Gl2i 



where If 
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where Af(x,y) 



7T 



v A* 



m 



f 



M 2 H xy 



f=u, d, s , c, t , b , e , jj, , r 



f 

1 fl-A 

dx 



dy 



1 — Axy 
A/ ( x, y) 



(3.1.6) 



(3.1.7) 



The integral If is evaluated in section A-3 of the appendix (see Eq. (|A-3-5| ) - Eq. ( A-3-13| ) ). It turns 
out that, due to the dependence of the factor for the H — > / / vertex coupling on the mass of the fermion 
nif, that the only fermion triangle which gives a significant contribution to the amplitude is for the case 
where mf 3> Mh- From the table, this is true only for the top quark mj. Hence, it turns out that the 
only fermion triangle that is necessary to take into, is the top quark triangle, and the contributions from 
the rest of the triangle sub-processes formed by the rest of the quarks, and the leptons can be neglected. 
Using this result, the amplitude for the contribution of the fermion triangle has the expression 



A 1 

Af = A f (« - 9^(qi • q 2 )) where Af = A G F 2 ~ 4a ™ (3 .1. 8 ) 
3.2 Boson loop sub-processes for Higgs production in 77 fusion 

For central exclusive Higgs production, the Higgs can also be produced through the subprocess 77 — > 
boson loop — > H, where the possible boson loops are shown in Fig. |B| ( taken from ref. Q). H~ in Fig. || 
is an un-physical, charged Higgs boson, and (f)^ is a Fadeev-Popov ghost. The formalism for calculating 
the amplitude of each diagram in Fig. [|, is similar to the approach used to calculate the amplitude of the 
fermion triangle above, in section [O]. Similarly here, after integration over the unknown momentum I in 
the loop of each diagram in Fig. || and after integration over Feynman parametres, the expression for each 
diagram takes the general form Q 
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O:, 



2tt 



2 



rat 



2m?„ 



(3.2.1) 



where d is the dimension of space-time. Terms proportional to and q\ were assumed to vanish, and 
from the kinematics shown in Fig. ||, it was assumed that q\ ■ q\ = ■ In the limit that d — > 4, the 
term proportional to T ( 2 — |) on the RHS of Eq. ( |3.2.1| ) tends to infinity. However, when one sums over 
the contributions to the amplitude given by all the boson sub-processes shown in Fig. ||, these divergencies 
cancel exactly (see table below). One requires also, that this sum over all the amplitudes for the boson 
loop sub-processes shown in Fig. ||, satisfies the condition 



(3.2.2) 



such that the amplitude for the sum is gauge invariant. In ref. f7|, this sum was taken and the result was 
an almost gauge invariant expression, since terms proportional to v# and higher were neglected. In the 
calculation which lead to the results in this paper, the result gives an exactly gauge invariant expression, 
after using dimensional regularisation to remove terms which do not satisfy the gauge invariance condition. 
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Hence, plugging the results shown in the table for the coefficients into Eq. ( 3.2.1| ), the result of taking 
the sum over the amplitudes for the sub-processes shown in Fig. || gives a gauge invariant expression, and 
in the limit that d — > 4, the divergencies cancel exactly, such that the expression of Eq. ( 3.2.1 ) reduces to 



A h 



v A 1 



where Ah 



a, 



7T 



(3.2.3) 
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Figure 5: Boson loop subprocess for Higgs production in 77 fusion 



It should be noted from the results for the amplitude of the subprocess of Fig. [| (a), the subprocess 
77 — > W triangle — > H, interferes destructively with the subprocess 77 — ► fermion triangle — ► H shown 
in Fig. f§. 
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3.3 The cross section for central exclusive Higgs production through 77 fusion 

Now that the amplitudes for the sub-processes 77 — > fermion triangle — > H and 77 — ► boson loop — * H 
have been calculated, the results can be plugged into Eq. ( fO| ) to derive the cross section, for central 
exclusive Higgs production through 77 fusion. The result, taking into account all possible sub-processes 
shown in Fig. || and Fig. [| is found to be 

<r^ c (p+p^ p+H+p) = 0.1 fb (3.3.1) 

4. Conclusion 

The results of this paper are summarized in the table below. <r exc is the exclusive cross section, which 
includes multiplication by a factor for the survival probability, for central exclusive Higgs production. The 
results are given for the mechanisms pp — > 77 — > H and pp — > PP — ► H. Note that in these results, the 
cross section for central exclusive Higgs production in the case of 77 fusion, is multiplied by a factor for the 
survival probability of 1. This is because in the case of photon exchange, there is no hard re-scattering to 
suppress, and the large rapidity gaps between the Higgs and the two emerging protons are automatically 
present. 



process 


< | S 2 | > 


(T exc (fb) 


PP 


0.023 


2.7 


PP 


0.004 


0.47 


77 


1 


0.1 



The results show that, taking the survival probability to be 0.02, which is the value used in ref. 0, then 
the result for crf^p for central exclusive Higgs production at the LHC, almost agrees with the prediction of 
ref. (which was 3fb). However, if the survival probability is an order of magnitude smaller as predicted 
in ref. Q, then cr*jj?p will be an order of magnitude smaller and, it becomes competitive with c^ xc for 
central exclusive Higgs production at the LHC. 
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A. Appendix 



A-l Evaluation of the integral over the anomalous dimensions 71 and 72 of the momentum 
Q 2 in the gluon density function 



The Born amplitude was calculated in Eq. ( |2.8| ) , in terms of the gluon density as a function of the anomalous 
dimensions 71 and 72, for the two gluon ladders in Fig. |l|. One now needs to integrate over 71 and 72, and 
also over the momentum in the t-channel gluon, namely Q. Altogether the necessary integrations take the 
form 

Mjpjp (p+p^ p+H+p) =,47r 3 S J d 7l d l2 J ^e- s ( k ^ E l)f(x u Q 2 )f{x 2 ,Q 2 ) (A-l-1) 

Firstly, the integral over Q± is evaluated using the steepest descent technique. The gluon density is 
given by 

/ (Q 2 ,^ 2 ) = 2 (Q 2 ) 71 ' 2 e^ ln ^T (A-l-2) 

where In ( ~ ha ( 1 where so ~ lGeV. This comes from the BFKL ladder gluon exchange 

(see Fig. |l]), while the coefficient was taken from MRST - NLO - 2002 data (see Ref/[?]). u> (ji^) is the 
BFKL kernel defined as 



w (71,2) = ol s x (71,2) = ol s (tp ( 1) - V ( 71,2) - i> ( 1 - 71,2) ) (A-l-3) 

where ip ( /) is the digamma function and V ( /) = ^p~- In Eq. ( |Ad^ ), S ( k 2 ± , E 2 ) is the Sudakov 

q factor with the typical value ||, ||| S (Q\, E 
Using this substitution Eq. ( |A-1-1 ) then becomes 



form factor with the typical value [|, £§ S (Q 2 ± ,E ± ) = ^ (in (^)) , in the notation that E ± = ^f. 



/CO p 
^71^72 / 
OO J 



dQ 2 



Ql 

xexp(-^(Qi))expfa;(7i) In ^ + w( 72 ) In ^\ (A-l-4) 
where <f> (Q\) = ^ (in (^Jj - ( 7l + 72 - 1) hiQi (A-l-5) 



Differentiating the right hand side of Eq. ( A-l-5 ) with respect to In Q 2 ,, one sees that (j) has a saddle 



point at In Q\ = In • 



^ + £(71 + 72 



1). Hence, changing the integration variable to u 



In Q 2 ^, and 



expanding (p around the saddle point, Eq. (|A-1-4| ) can be written as 
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d 7l d 72 / dueV * {U ~ U0) ~^~) er {7l) ln ^^(t^) 

where ii = In + J^- (71 + 72 - 1) (A- 1-7) 

4 3a, 



Now the right hand side of Eq. ( A-l-6| ) has reduced to a Gaussian integral over u, which can be 



evaluated by the steepest descent technique, to give the expression 

Mjpjp ( p+p -> p+H+p) =AAir A s (^-) * J ^72 exp ^( 7l + 72 - 1) ( 7l + 72 - 1) + ln \ 

x exp ( uj ( 71) ln — + a; ( 72) ln — ] (A-l-8 

The BFKL function w (7) has a saddle point at 7 = |. Near to this point u> (7) can be written as 

1\ 1 / 1\ 2 .. /I 



^(7i, 2 )=^2 1 + 2(71,2-2) (A-l-9) 



Hence using Eq. ( A-l-9| ), Eq. ( [A-l-8 ) can be reduced to 



Mj>j> ( p+p -> p+H+p) =4^vr 4 s f ^-^ 2 y d 7 ic?72 exp ( / ( 71, 72) ) (A- 1-10) 
where the function / ( 71, 72) has the form 

/(71 , 72) ^(^f-^( 71 _^V(i) 1 n^^( 72 -i)V(i) 1 ,fi 

+ (71 +72-1) (jL ( 71 + 72 - 1) + ln ^\ (A-l-11) 
This function has a saddle point with respect to 71 given by 

7r= -fe(*-i)--ff + K(i)-S (AJJ2) 

(fe+^'fa^s) 



Hence, expanding / ( 71, 72) around 7 * p , the integration over 71 is evaluated using the steepest descent 
technique to give the expression 
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■ sp 



M PP ( p +p - p +H+P ) =2A,'., (£-)*[ d 72 , (A-l-13) 

" fe+1 »i? +u .(i)i» s : 



Now the function / (71, 72) has a saddle point with respect to 72 given by 



^2 



1 M H 



In 5i 

so 



2?r 
3a s 



In 



+ W " (I) In 5i 



(A-l-14) 



Here in the second line it is assumed that s\ ~ s^. For large s\ and S2, * s approximately ^. Using 
the same method as above, expanding /( 7^,72) around 7^, the integration over 72 is evaluated using 



the steepest descent technique for Eq. (A-l-13), to give the result 



Mjpjp ( p+p -> p+H+p) 



where 



/(7f,7f) 




-2exp( / ( 7 f,72 P -^)) 

M~ 



^+ln^+o;"(i)ln^ 



3a s 



a; 



(A-l-15) 



(A-l-16) 
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A-2 Feynman rules for the standard electro-weak theory 




H(q) 



Figure 6: 



\ 9 Y9» 

=-(i/2)ig,M„ \ =(+/-)iep n 

yw w .-WWW H » 



+/- 



Figure 7: Feynman rules in the standard electroweak theory 

A-3 Evaluation of the integral over the momentum in the fermion triangle loop 

The amplitude for the fermion triangle subprocess, summed over all quark flavours q = (u , d , s , c , t , b 
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and lepton flavours L = ( e , fi , r ) , for the mechanism 77 fermion triangle — > H was found in equation 



Eq. (3.1.4) to take the form 

i r d d l f 1 [ l ~ x 8m)((q^-^g^)(l-2xy)-2xy^fg^+rn)g^ 



(2tt) Jo Jo (p-A^x^) 



3 



+ 2y"4vrQ em (G F \/2) 2 / - — / cte / dy—^ where 

(2tt) Jo Jo (fa - A/ (*,!/) 



/ 

A/ = mJ-M&xy and ^ = ^ + ^ (A-3-1) 

/ <J L = e,fj,,r 

where ^2 q denotes the sum over all six quark flavours q = u , d , s , c , t , b and ^2 L = e „ T denotes the 
sum over all three lepton flavours L = e , \x , r. 



From the numerator of the integrand, on the RHS of Eq. ( A-3- 1| ), one can see that there is a gauge 

t ■■ M 2 \ - 



invariant term yq^q^ ^9^") (1 — 4scy) , and the numerator in the integrand of the second line on 

the RHS gives a vanishing contribution to the integration over I, for d — > 4. However one is still left 
with the terms —2xy—j L g iiV + m 2 g^ v in the numerator of the integrand, which are certainly not gauge 
invariant. However, conveniently this non gauge invariant piece is exactly equal to Ay (x,y) introduced 
in Eq. ( 3.1.2; ). To deal with this non gauge invariant piece, it is useful to use dimensional regularisation, 



when integrating over the I 2 term in the numerator of the integrand on the second line. In this approach, 
one initially integrates over d + e dimensions in the limit that e — ► and d — > 4. In this way the non gauge 



invariant terms disappear. Hence, evaluating the integral over I on the RHS of Eq. ( |A-3-l| ) gives 



A>t u = 2 V4vra em G^24^ > \ *L / dx / dy^±± j— , 

1 f r(3)(4vr)2 Jo Jo &f(x,y) 

x 1 ( -\) 3+1 Y (3 - 4+£ - l) /•! fl-x 

+ 2 lim y 4vrQ em G|23 i '- K - §- L dx dy 8m) ( AgT ~(d + e) O (A-3-2 

£ -*°f ' 2r(3)(4^)— Jo Jo 



In the limit that e — > and d — ► 4, the gamma function T ( 3 — — l) — ► ~ and the RHS of 
Eq. ( |A-3-2| ) reduces to 

i 1 1 ,1 f i-x 8m) ( ( - *%.gT) (l — 4xy) +Af(x,y) gT) 

1 j 2(4vr) 2 J Jo A f (x,y) 

+ 2 V4vra em G}23-— / dx f dy8m)eg^ (A-3-3) 
, e2(47r) Jo Jo 
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Thus, after canceling e in the numerator and the denominator in the second line on the RHS of Eq. (A 



3-3 |), th e second line exactly cancels the non gauge invariant part of the integrand on the first line. Hence, 



one is left with the purely gauge invariant expression 



K = - 2 ( « - ^r) E *f (A-3-4) 



where If is the only remaining integral to evaluate, which takes the form 

E If = E mj J dx J dy ^ ^ V where A f (x,y) = mj - M 2 H xy 

f f=u,d,s,c,t,b,e,fji,T ^ ' 

(A-3-5) 

To evaluate this integral, this are two cases to consider, namely (1) when mj 3> Mj^, which is true 
for the top quark when m/ = mt = 175 GeV, and (2) when mj <C M|, which is true for all the rest of 
the fermions listed in the table. Therefore Y^f If can be separated into two parts, namely 



E J / = J * + E J / (A-3-6) 
For case (1), where mj = mt S> Mh , one can write It in a more convenient way as 



m? f 1 , f l - x , 1 - 4xy m 2 t f 1 , f 1 ^ , / 1 " 4 m 

It = T7T I dx I = T7T dx d V 4 + ~ 

, . ii .ii v T, ;( 



M' H Jo Jo Ii£ - X y M H Jo Jo \ - X y 



mi f ( „ A m 2 \ f 1 dx , / Af£ , 

i" 4 ^ / "In l--fx(l-x) 



M 2 H V V io z V m? 



^ ' 2 + l-4 m 



€Hi (i -*> + Kf) 2 * (i -* )2+ Ki) 3 ^-^ )) 



(A-3-7) 



where in the last step, the logarithm was expanded in a Taylor series around x = 0. Evaluating the 
integral over x, and since it is assumed that M\\ <C tti^ , retaining terms no smaller than ( — ) , the RHS 



of Eq. ( A-3-7| ) becomes 



I (l+JL** + ....) (A-3-8) 
3 l 120 m ' 
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For case (2), where m/ <C Mh which includes all the fermions in the table except for the top quark, 
there are two possible regions of integration, namely (I) when Ml x y > m 2 and (II) when Ml xy < m 2 . 



In the region where Mj^xy > mj , the RHS of Eq. ( A-3-5| ) reduces to 



^2 jregion (I) 



m^ 
^ Ml 



E 



f 

Ml 



H 



,2 dx 
f 

1 



1-x 



dy 



1 — Axy 

xy 



dx , 
,2 — hi 

f x 



M 



x(l 



X 



m 



f 



5^ 9 /if 2 



m 2 j l 



m 2 f . m 2 f ( 

+ TTTPO^og ( 2 , x = 1) - 2^ V^-PO^og 2 , ; 



mj Mfj f£M B 



Ml 



+2 y i + 4 

^ 1 Ml 



m. 



Ml 



In 



m. 



Ml 



(A-3-9) 



In the region (II) where M^xy < mi , the RHS of Eq. (A-3-5) reduces to 



E region (II) 



( 1 - 4xy) 




E 



% Mix 



'Ml 



(A-3-10) 



Hence, adding the contributions of Eq. ( A-3-9|) and Eq. ( A-3-10| ) for the contributions of region (I) and 
region (II) of the integral, gives the result 



E'/ 



Ml 

2 

m, 



1 - 6 



f 

Ml 



1. Ml 



lu- 



ll 



m, 



In— ^ 

mi 



[ polylog ( 2 , rc = 1) - polylog I 2 , rc = 



/?? 



for mj < M// 



/ 



+ 2 



(A-3-11) 
(A-3-12) 



From Eq. ( A-3-12 ) and Eq. ( A-3-q ), the result for the evaluation of the integral If has its main 
contribution from the top quark triangle, such that 
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If ~ It , = ~ for mt » Mjy (A-3-13) 
/ 3 

Plugging this result into Eq. (A-3-4) gives the final expression for the amplitude of the fermion triangle 
subprocess shown in Fig. ||, for the sum over all quark q = (u , d , s , c , t , b) contributions and lepton 
contributions L = e , \i , r as 



Af = A f U$ - ^V") where A, = - 2 -^H^l {A . 3 . U) 
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